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Experimental results and their interpretations are presented on the nonlinear acoustic effects of 
multiple scattered elastic waves in unconsolidated granular media. Short wave packets with a central 
frequency higher than the so-called cut-off frequency of the medium are emitted at one side of the 
statically stressed slab of glass beads and received at the other side after multiple scattering and 
nonlinear interactions. Typical signals are strongly distorted compared to their initially radiated 
shape both due to nonlinearity and scattering. It is shown that acoustic waves with a deformation 
amplitude much lower than the mean static deformation of the contacts in the medium can modify 
the elastic properties of the medium, especially for the weak contact skeleton part. This addresses 
the problem of reproducibility of granular structures during and after acoustic excitation, which is 
necessary to understand in the non destructive testing of the elastic properties of granular media 
by acoustic methods. Coda signal analysis is shown to be a powerful time-resolved tool to monitor 
slight modifications in the elastic response of an unconsolidated granular structure. 

PACS numbers: Valid PACS appear here 



I. INTRODUCTION 



Unconsolidated granular materials are known to ex- 
hibit a high complexity in their macroscopic behaviors. 
This gives rise to numerous fundamental and applied 
processes that are currently intensively studied, such as 
avalanches, dune formation, compaction etc [l|, [2, H, Hj] . 
Most of the studied physical processes concerning the un- 
consolidated granular matter are obviously related to the 
possibility for the grains to move relative to each other, 
from microscopic (much smaller than a grain diameter) 
to macroscopic (much larger than a grain diameter) rear- 
rangements. While large movements of grains have been 
observed and modelled for the processes where the gran- 
ular material behavior may be fluid-like (even complex 
and not classical), like avalanches, microscopic relative 
movements are less intuitive and less studied. They are 
however the only possible relative movements in solid-like 
granular media where a static stress is applied and the 
grains are confined in a given volume. 

In this context, the acoustic waves, known to be highly 
sensitive to the contact stresses 0, H S 13 : can be useful 
for the monitoring of these small relative movements Q , 
impossible to detect by other existing methods. They 
could even be useful for the generation of these small 
movements. Among the several existing methods for 
the grain movements or contact stresses monitoring (car- 
bon paper experiment, photo-elasticity, image processing 
from CCD camera, X-ray tomography or photo-imaging 
in iso-index configuration 

[E Sf, E3, G3]), the acoustic 
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methods are, to the best of our knowledge, the only ones 
sensitive to the contact stresses inside the volume of a 
3-D non transparent medium. One of the main advan- 
tage of acoustic methods is also the temporal resolution 
of the measurement, which could decrease down to the 
wave period (one millisecond at one kHz), and is partic- 
ularly interesting in transient processes like avalanches 
for instance. However, there is currently a lack of un- 
derstanding of the acoustic transport properties through 
granular media, which makes these methods mostly qual- 
itative. 

Several previous works @, [H, [H[ have shown the 
following qualitative features for the acoustic wave prop- 
agation through stressed and non cohesive granular me- 
dia. The propagation of a low amplitude (with a strain 
much lower than the average static strain of the con- 
tacts) and low frequency (with a wavelength much larger 
than the grain diameter) wave can be linear and is mostly 
governed by the averaged properties of the medium (aver- 
age density, average static stress, coordination number) 
[3; EH ■ The propagation of a high amplitude (with a 
strain lower than the average static strain of the contacts 
but which can be comparable) and low frequency wave is 
nonlinear and governed partly by the so-called weak con- 
tacts (with a static strain lower than the average one) 
0, EI 113 • This is a consequence of the fact that in- 
dividual stress-strain relationships for each contact have 
higher values of nonlinear parameters when the pre-stress 
is smaller. For acoustic waves with wavelengths compara- 
ble to the bead diameter, however, these simple features 
are not straightforward, and the transport properties de- 
termination remains a fundamental problem. How to de- 
scribe the wave scattering in a network of beads with 
geometrical and contact disorder? What is the role of 
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the so-called force chains in the transport properties of 
short wavelength acoustic waves? Can the weakest con- 
tacts make the multiple scattering process nonlinear? 

A linear diffusion approach has been recently applied 
to the analysis of coda-type signals in confined granular 
media ■ This is one of the first attempt to understand 
the transport properties of short wavelength elastic waves 
in such media at the laboratory scale and to extract pa- 
rameters such as a characteristic time of absorption r a 
and a diffusivity D. This approach is based on previous 
works in elastic wave diffusion carried out for slightly 
different media (solid rods in water [lj|, chaotic cavities 
or glass bead in water [2(| for instance) and on a 
different scale pll^. 

In this article, we report some experimental results of 
nonlinear acoustic wave propagation, with initially radi- 
ated wavelengths close to the bead diameter. In some 
well-chosen experimental configurations, it is possible to 
observe coda-type signals, together with a nonlinear sig- 
nal of low frequency, which is identified here as a self- 
demodulated signal. For some widely encountered config- 
urations, linear and nonlinear wave processes of multiple 
scattering, self-demodulation and self-action are observed 
[23l |24| . The last nonlinear process of self-action of the 
multiple scattered field is observed even for moderate ex- 
citation amplitudes (deformations much lower than the 
average static strain of the contacts). The application of 
some of the presented effects to the characterization of 
granular media is finally discussed. 



II. EXPERIMENTS 

As a preamble, carrying quantitative experiments on 
acoustic wave propagation in unconsolidated granular 
media is a hard task for several reasons that are described 
in the following. 

First, in order to excite the modes of propagation 
inside the solid frame, i.e. through the elastic beads 
and their contacts, and not through the fluid saturat- 
ing phase, it is necessary to put the piezo-transducer in 
direct contact with the medium itself (or in contact with 
a solid which is itself in direct contact with the medium) . 
Other types of transducers for the excitation of acoustic 
waves, like laser generation of ultrasound for instance, 
have not been applied to our knowledge, certainly be- 
cause of the high efficiency needed in these strongly at- 
tenuating media. The application of a mechanical stress 
with a contact transducer has several consequences on 
the possible experiments. What can be done usually in 
liquids for acoustic metrology, i.e. changing the emitter- 
receiver distance in order to measure the wave velocity 
for instance, is extremely difficult in unconsolidated gran- 
ular media, because the contact network is modified by 
any variation of external conditions. 

The second difficulty in carrying quantitative experi- 
ments is the experimentally observed strong influence of 
the applied static stress on the transmission coefficient es- 
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FIG. 1: Schematics of the experimental configuration. 



pecially for short wavelength acoustic waves. This result 
is discussed in the following. The consequence is the im- 
possibility to perform quantitative reproducible and pre- 
cise experiments if ones modify the applied static stress 
on the sample, even if it is measured to be equal after 
a back and forth modification. Moreover, if one would 
be able to apply exactly the same stress, the memory 
effects due for instance to the hysteresis in the quasi- 
static stress-strain relationship, take place [25|, |2f3, [23. In 
an attempt to control carefully the granular preparation, 
the acoustics of compaction has been recently studied 
[28| . However, it allows to perform experiments only for 
weakly stressed granular packings (typically from 75 Pa 
to 3 kPa), far from the stress necessary for the observa- 
tion of the signals described here. 

The third difficulty, is related to the fragility of such 
media [H, |j[ ■ Due to the inhomogeneity of the individ- 
ual contact loads EH EH , some of the contacts 
are weakly loaded, which make them able to clap, slip or 
slide under the action of a small amplitude strain wave. 
In contrast, the so-called strong contacts, are supposed to 
take part in the force chains which constitute the elastic 
skeleton of the bead arrangement. Previous works show 
by simple analytical formulae or experimentally that the 
weaker the contact the higher is its acoustic nonlinearity 
Q. It is thus expected that the major role in the lin- 
ear acoustic properties is played by the strong contacts 
and the force chains while the nonlinear acoustic prop- 
erties are more controlled by the weak contacts. There 
exist few experimental configurations (in particular when 
short wavelengths are used) when these roles cannot be 
qualitatively separated like that, and either weak con- 
tacts influence importantly the linear properties [3fj| or 
strong contacts are involved in the observed nonlinear 
effects. 
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A. Experimental setup 



In Fig. [H a schematic representation of the experimen- 
tal setup is shown. A wave packet made of a sine wave at 
100 kHz modulated in amplitude by a Gaussian function 
is high-pass filtered over 40 kHz (in order to avoid any di- 
rect low- frequency excitation) , amplified and launched in 
the medium with a shear piezo-transducer. Due to mul- 
tiple scattering, mode conversions (including dilatancy), 
the acoustic energy is rapidly distributed between longi- 
tudinal and shear vibrations in the medium. The use of 
a shear emitter has the advantage of avoiding the direct 
excitation of a pressure wave in the air saturating the 
beads. In the following results, two types of Gaussian 
function widths for the initially launched wave packets 
have been used, and two types of surface treatment for 
the shear transducer (the first type is the raw transducer 
membrane with asperities estimated to be of character- 
istic scale ~ 50 /jra and the second type is a transducer 
with a glued single layer of 2 mm diameter glass beads, 
the same beads as in the medium itself). The container 
diameter is D = 20 cm and the propagation length ranges 
between 1 cm and 18 cm. The applied static uniaxial 
stress is measured with a force sensor at one end of the 
container and ranges from 30 kPa to 500 kPa. A longi- 
tudinal piezo-electric receiver is placed at the bottom of 
the container. This transducer has been chosen because 
of its high sensitivity over a wide-frequency band, allow- 
ing for a large range in the acoustic amplitude measure- 
ments. In all the experiments presented in the following, 
the estimated average contact strain Eq (of the typical 
range 3 — 7 x 10 -4 for a 200 kPa applied stress) is always 
greater than the maximum acoustic strain e max ~ 10~ 5 . 

When using this setup, it is possible to modify the 
emitted acoustic frequency, the acoustic amplitude, the 
applied static stress, the propagation distance. But im- 
portantly, it is better for a quantitative insight in the 
wave propagation phenomena, to modify parameters of 
the setup that change the smallest number of medium 
(or wave) parameters. For instance, changing the fre- 
quency or the static stress modifies the wave scattering, 
absorption, dispersion and nonlincarity. Changing the 
distance, changes irreversibly and in an uncontrolled way 
the granular medium state, as it is necessary to unload 
completely the medium and add or remove some granu- 
lar matter. The convenient parameter to modify is the 
wave excitation amplitude, which can be accurately mea- 
sured, and is reproducible. The medium is a priori not 
irreversibly modified. 

In the following several experimental results are pre- 
sented and their implications on the granular material 
elasticity is discussed. The parameters of the experiment 
may be different, they are consequently recalled for each 
group of results when necessary. 
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FIG. 2: Experimental temporal signals that are typically re- 
ceived for different excitation amplitudes. The applied static 
stress is ~ 300 kPa and the emitter - receiver distance is 
D = 10.5 cm. 



B. Multiple scattering and nonlinearity 

The signals presented in Fig. [5] are typical experimen- 
tal signals that can be observed when the initial acoustic 
wavelength is of the order or less than the bead size and 
for a static stress larger than ~ 100 kPa). For an estima- 
tion of the wavelength, one can use the elastic parameters 
of glass and the stress-strain relationship of an average 
loaded contact in an effective medium theory. However, 
it is important to keep in mind that this definition ne- 
glects wave velocity dispersion which can be important 
when the spatial scale of the beads coincides with the 
wavelength. As an estimation, the wave velocity in a dis- 
ordered three-dimensional packing of 2 mm glass beads 
with an applied static stress of 200 kPa is ~ 300±40 m/s, 
which gives a wavelength of ~ 3 ± 0.4 mm at 100 kHz. 
This means that for a wave packet centered on 100 kHz, 
submitted to a 200 kPa static stress, the wavelength is of 
the order of the bead diameter, and consequently, strong 
scattering occurs. The typical experimental signals re- 
ceived in this case are strongly distorted due to multiple 
scattering, lasting from 10 to 100 times longer than the 
emitted signal (see Fig. This feature is widely ob- 
served in seismics [2l1 | or in multiple scattering experi- 
ments, and has been observed in roughly the same con- 
ditions in glass bead assemblage [l§|. However, we found 
another feature in this case, not mentioned previously 
in the literature: for different excitation amplitudes, the 
normalized signals in Fig. [2] are different. This relative 
amplitude and shape dependence is the signature of a 
nonlinear process. One obvious difference between the 
traces is the emergence, for an increasing excitation am- 
plitude, of a signal at earlier times, relatively lower in 
frequency than the 100 kHz coda-type signal. We veri- 
fied in several other experiments 0, [HI that this LF sig- 
nal is the self-demodulated contribution associated with 
the rectification (demodulation) of the initial HF wave 
packet. Its nonlinear nature is evidenced by its nonlinear 
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FIG. 3: Experimental spectra corresponding to temporal sig- 
nals of Fig. [2] The spectrum of the emitted signal with a 22% 
FWHM (full width at half maximum) is superimposed on the 
top curve. 



amplitude dependence on the excitation amplitude. This 
is an important difference compared to results [H, H3 |. 
where the very same types of signals were observed (for 
instance for e a = 5 x 1(T 6 in Fig. [2]). In [3, the LF 
contribution is identified as the linear coherent part of 
the propagated initial pulse, and the HF contribution is 
interpreted as the multiply scattered signal part. A pos- 
sible explanation for such a different observation (mainly 
in the amplitude dynamics behavior and physical nature 
of the LF signal) is the difference in the spectral width 
of the emitted signal spectra. In Ref. [12||, the signal 
spectrum width is relatively wide, which may allow for 
the direct radiation of the linear coherent LF component 
in the medium. Then, the observation may be inter- 
preted as a strong frequency dependent transmission in 
the medium (due to absorption and scattering), which is 
able to modify importantly the transport character of dif- 
ferent frequencies (propagative as in an effective medium, 
multiply scattered, diffusive...). 

In the present case however, when the emitted spec- 
trum is sharper, the direct radiation of the LF signal 
is negligible (and in any case is avoided by a high-pass 
filter in the setup). This LF signal is generated in the 
medium itself, by nonlinear processes of frequency mixing 
, widely described in the literature [HI, [H, [H, [H[ . 
This can be seen in Fig. [31 where the relative amplitudes 
of the LF and the HF contributions are modified by in- 
creasing the excitation amplitude. 

In Fig. [31 the received temporal signals are presented 
for three source-receiver distances, with the same ap- 
plied static stress 300 kPa (measured with a precision of 
~ 3 % which is not sufficient to perform precise quanti- 
tative measurements). An important qualitative feature 
observed on these normalized (by their maximum) traces 
is the relative increase of the LF contribution compared 
to the HF one with distance. This feature is usually ob- 
served for the nonlinear self-demodulation process [3l| . 
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FIG. 4: Signals received after propagation through different 
thicknesses of granular medium, for identical static pressure 
of 300 kPa and excitation amplitude £ a =5xl0~ 6 . 



where over a sufficient distance, only the LF contribution 
is detectable, due to the difference in attenuation of the 
HF and the LF waves. In this experiment, the LF signal 
absolute energy is measured to be roughly constant with 
distance, while the received HF signal energy is decreas- 
ing drastically. The former is a result of the competition 
between nonlinear effects, which pump low frequencies 
through frequency-down-conversion processes, and linear 
attenuation. 

In Fig. [SI the total energy of different signals is plotted 
as a function of the applied static stress for a propaga- 
tion distance D = 12.5 cm and an excitation amplitude 
£ a — £max/50, where Emax — 10~ 5 . Three signal energies 
are defined here: the energy of the HF contribution Ehf, 
the energy of the LF contribution E^p and the energy 
of a directly radiated LF signal E l LF (in this case, a LF 
pulse with the same frequency content as the nonlinearly 
demodulated pulse is directly excited by the emitter and 
received after linear propagation). First, comparing the 
received energies Ehf and E l LF evolution as a function 
of the applied static stress, an important slope difference 
is clearly visible: E l LF dependence on stress is less than 
to a power 1 and Ehf dependence is close to a power 4 
of the stress. In order to explain this difference in the 
evolution of the transmitted energy as a function of the 
applied static stress, one should recall that two dominant 
processes of sound attenuation may play a role: the scat- 
tering and the absorption due to linear processes. Con- 
sidering that the scattering of the HF ~ 100 kHz, is much 
more important than for the LF ~ 10 kHz, the observed 
difference may be attributed to the strong stress depen- 
dent scattering. This is consistent with the simple idea 
that by increasing the static stress, contacts are created 
which adds new paths (or strengthen existing paths) of 
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FIG. 5: Signal energy dynamics as a function of the applied 
static pressure. The HF (Ehf) and LF (E2 F ) contributions 
are separated by filtering. The energy of a reference LF (E l LF ) 
linear signal is also plotted for comparison. 



acoustic energy transmission. Moreover, waves follow- 
ing force chains with a given stress may change from 
an evanescent character to a propagative one with in- 
creasing stress, because the cut-off frequency (sometimes 
called the Einstein frequency) is increased [3(J Hfl H3] ■ It 
is then possible to ask if this scattering which strongly 
depends on static stress depends also on the acoustic am- 
plitude of excitation. Is the acoustic wave able to switch 
dynamically propagation paths or force chains? The aim 
of the next section is to analyze this opportunity and its 
possible manifestations. 

The competition between nonlinear effects and atten- 
uation is also visible in Fig. [5] when comparing the dy- 
namics of E e LF and E L l F . For signals having the same 
low frequency propagating through the same medium but 
generated cither directly by a transducer or by nonlinear 
effects inside the medium itself, the energy dynamics as 
a function of the applied static stress is different. It can 
be assumed that the linear attenuation varies in the same 
way for both signals while the nonlinearity of the medium 
and the HF attenuation both playing a role in the genera- 
tion of the LF signal 3l| may be strongly affected. These 
modifications only affect the nonlinearly generated sig- 
nal. The E LF dynamics as a function of the applied static 
stress is almost constant, while the E l LF is increasing. It 
means that the efficiency of the nonlinear interactions 
leading to the self-demodulated signal is decreasing with 
increasing static stress. Two main processes can lead to 
a diminishing of the nonlinear interaction efficiency: the 
HF attenuation increase (the amplitude of the nonlinear 
sources diminishes) and the intrinsic nonlinearity of the 
medium decrease. As the HF attenuation (probed by the 
dynamics of Ehf) is strongly diminishing with increas- 
ing static stress, it means that the nonlinearity of the 



medium is decreasing with increasing static stress. This 
is consistent with estimates based on the Hertzian con- 
tact nonlinearity and experimental observations @, EH . 

Now, we analyze the temporal features of the HF sig- 
nals, namely the codas. Coda features can be analyzed 
by performing an envelope detection (for instance by in- 
tegrating the square field over a sliding window). This 
may cause problems for quantitative analysis because one 
has a degree of freedom for the size of the sliding win- 
dow. In principle, the diffusion approximation is shown 
to be valid for the configurational average intensity (/) 
of the incoherent field under some conditions [38|. The 
incoherent field can be defined as the contribution which 
cancels out with averaging over different realizations of 
the disorder in the medium [H, H^] . The average inten- 
sity (/) is usually accessed experimentally by the use of 
an intensity sensor, which is large compared to the wave- 
length (a photo- multiplier in optics for instance). For 
elastic waves however, the piezo-electric transducers are 
sensitive to the field itself. It is thus necessary, in order to 
record (I), to record wave field signals, then to calculate 
the corresponding intensities, and finally to average the 
different obtained intensities over different configurations 
of the disorder. However, when the receiving transducer 
has a sensitive surface with dimensions larger than the 
wavelength, a spatial averaging is performed on the field. 
In this operation, the incoherent contribution is dimin- 
ished compared to the coherent one. Consequently, even 
if most of the energy is incoherent at the receiver loca- 
tion, there could be some situations where the detected 
signal is composed of coherent and incoherent parts of 
comparable importance. 

In unconsolidated granular media, the configurational 
averaging can be performed by several means, which may 
not be equivalent. One could use arrays of transduc- 
ers for the excitation and detection of the multiple scat- 
tered waves, each transducer being small compared to the 
wavelength. Then without modifying the medium, each 
emitter - receiver couple, for the same distance through 
the medium could provide a signal to be treated. Other- 
wise, with one emitter and one receiver, one could mod- 
ify the structure of the medium (geometrical and contact 
force disorder) by keeping constant as much as possi- 
ble the macroscopic external conditions @. These two 
means can allow in some limits for an estimation of the 
diffusivity and the absorption characteristic time of the 
diffusion equation. The first method necessitates small 
enough transducers to avoid averaging of the field, which 
causes problems of coupling with the medium when the 
wavelength is of the order of the bead diameter and the 
transducer size. Some of the transducers may not be in 
efficient contact with beads. The second method, due to 
the memory effects [U SO] or the drastic influence of the 
external parameters on the elastic properties (see Fig. [5] 
for instance), does not guarantee to be in the same ex- 
perimental conditions from one realization of the disorder 
to another. Small perturbations or local modifications of 
the contacts forces in the medium are consequently diffi- 
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FIG. 6: Typical temporal signals (normalized by the excita- 
tion amplitude) that are received at a distance of 6 cm from 
the emitter for different excitation amplitudes. The emitted 
signal is a Gaussian wave packet centered on 100 kHz with a 
full width at half maximum (FWHM) of 6% in frequency, a) 
Row signals, b) High-pass filtered signals. 



cult to probe using averaged quantities, obviously more 
adequate for an estimation of average properties of the 
medium. 



C. Analysis of the amplitude dependent codas 

In this section, we investigate fast dynamic modifica- 
tions of the medium submitted to an acoustic wave, as 
well as slow dynamic modifications or permanent mod- 
ifications of the medium. The experimental configura- 
tion is slightly different from the one in the previous sec- 
tion. Some glass beads have been glued to the shear 
emitting transducer surface, in order to couple better to 
the medium and to have a better signal to noise ratio. 



The analysis of the obtained spectra indicates that it 
is impossible to detect neither the second nor the third 
harmonic component, certainly due to the strong ob- 
served attenuation above 100 kHz. However, the self- 
demodulated wave is quite easily generated and detected 
because of its lower attenuation. In the following, the HF 
coda contribution and this LF self-demodulated contri- 
bution are separated by post-processing filtering in order 
to analyze in details the different amplitude dependent 
processes. 

In Fig.^a), typical signals that are received by a sensi- 
tive wide-band longitudinal transducer (with a sensitive 
surface of 4.5 cm in diameter) are presented. The ex- 
perimental protocol elaborated in order to analyze the 
amplitude dependent effects of acoustic wave transport 
is presented in Fig. [TJd) . It contains 59 successive excita- 
tion amplitudes. The 10 first amplitudes, as well as the 
last 10 amplitudes, are identical and relatively low (the 
acoustic strain is e a = £1-10 — 10 -7 ). The odd num- 
bered amplitudes 11-49 are gradually increasing up to a 
maximum amplitude of e a ~ 10~ 5 , while the even ampli- 
tudes 10-50 are equal to the minimum probe amplitude 
e a ~ 10~ 7 . This allows to compare the response of the 
medium for different excitation amplitudes, but also, to 
visualize the medium modifications after being excited 
by a strong wave by using the weak probe wave. The 
experiment time for the acquisition of one signal of this 
protocol is close to 30 s. In Fig. [6l signals are numbered 
according to their excitation amplitude in the protocol, 
the first signal of the protocol being denoted by si. The 
excitation signal is a sine wave modulated in amplitude 
with a Gaussian function of 6% full width at half max- 
imum (FWHM) in the frequency domain, which, at the 
central frequency /o = 100 kHz gives a ~ 0.2 ms dura- 
tion pulse at half the maximum. This duration is roughly 
the time necessary for the wave to travel from the emitter 
to the receiver for the 6 cm distance. 

In Fig.[6{a), received signals over the full available fre- 
quency band (from 1 kHz to 300 kHz) are presented 
for several excitation amplitudes of the protocole: £10, 
en, £19, £20- The signal structure is quite complicated 
and shows several packets typically 4-5, having a shape 
close to the emitted signal. It is important to note here 
that the time dependent intensity of these signals is much 
different from a diffusive type enveloppe, even if the av- 
erage wavelength is of the order of the bead radius. An- 
other important feature is the emergence, at a sufficiently 
high excitation amplitude of a relatively low frequency 
(LF) signal (see for instance signal S19 corresponding 
to the amplitude £19 in Fig. [SJa)). Due to its ampli- 
tude dependent character (this signal is not visible for 
lower excitation amplitudes), this LF signal is nonlinear, 
and has been identified in previous works as being a self- 
demodulated contribution @, HH . 

In the following, we analyze in details the high fre- 
quency part of the signals, and for quantitative compar- 
ison, the raw signals are high-pass filtered (Fig. [HL)). 
These coda-type signals have comparable time durations 



to those of Fig. [2j but their envelope characteristic fre- 
quency is lower. This is a consequence of the narrower 
Gaussian spectrum of excitation. However, despite of 
the difference in coda shape, the characteristics of wave 
transport arc similar. The most important characteris- 
tics of the signal evolution with excitation amplitude is 
the fact that their temporal shape is strongly modified, 
even for very moderate changes in the excitation ampli- 
tude. This can be seen by comparing the shapes sio and 
sn in Fig. [DJb). Of course, by increasing further the 
excitation amplitude, the distortion is further increased 
(signal sig). When the excitation amplitude is decreased 
down to the lowest amplitude of the protocole, i.e. to 
the probe amplitude e\, the shape modification of S20 
compared to si is barely noticeable. In order to monitor 
quantitatively the modifications of the signal shapes due 
to amplitude dependent effects, we make use of the fol- 
lowing parameter R m , n defined from the inter-correlation 

function C n ,m{i) = Jq s n (t)s m (t — r)dr, where T is the 
observation window duration, by 



Rn.r. 



C n . m {0) 



(0) 



(i) 



The energy of signal s n (respectively s rn ) is propor- 
tional to C nn (0) (respectively C mm (0)). The normaliza- 
tion by y C nn (0)C mm (0) makes consequently the param- 
eter Rn.m independent of the actual energies of s n and 
s m , which is important in our case where signals with dif- 
ferent amplitudes are compared. The parameter is only 
dependent on the relative shapes of the signals and can be 
seen as the level of resemblance of two signals, equal to 1 
if the two signals are exactly the same and progressively 
diminishing when the signals differ. 

In Fig.[7]Ja), the parameter i? nj i is plotted for each sig- 
nal Si to Ssg of the experimental protocol. It is clearly 
visible that the parameter goes under a value of 

0.9 while the signals S2 to sio exhibit a very accurate 
correlation with s\ (R2— 10,1 — !)■ When the amplitude 
is increased further, the parameter R n ^\ decreases dras- 
tically to or less, which shows that a strong field mod- 
ification occured (the negative values of the parameter 
Rm,n do not have a particular meaning, except when the 
value reaches —1, which means that the signal phases are 
inverted). Concerning the even signals generated with a 
low excitation amplitude, their resemblance parameter 
remains practically equals to 1 up to the 30th signal of 
the protocol. This means that no modification of the 
granular medium elasticity is detectable by this acoustic 
probing. In contrary, when higher acoustic amplitudes 
are generated in the medium, the level of resemblance of 
the probe wave with the first signals of the protocol be- 
gins to fall down progressively, until the typical value of 
0.6 after the maximum excitation strain £49. The acous- 
tical properties of the medium have been modified by 
the strong wave action. A possible interpretation of this 
experimental observation could be found in the temper- 
ature effects on the linear and nonlinear hystcretic prop- 
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FIG. 7: (a) Resemblance parameter -R n> i between the signal 
number n and the signal number 1 of the protocol (d). (b) 
Resemblance parameter i? n>n _i between the signal number n 
and the signal number n — 1 of the protocol (d). (c) Arrival 
time T e of half the signal energy as a function of the signal 
number in the protocol (d). (d) Protocol for the successive 
excitation amplitudes of the experiment. 



erties of the medium [4l| . The wave action on the con- 
tacts increases the average temperature and consequently 
modifies the linear and nonlinear contact properties for 
some characteristic time which can be of the order of sev- 
eral acoustic periods. Elastic memory effects at the level 
of the contacts themselves (stick-slip motion for Hertz- 
Mindlin contacts for instance [I^]) or at the mesoscopic 
scale [13] could also contribute. 

When, at the end of the protocol, only the probe am- 
plitude is generated (signals S50 to S59), the parameter 
Rn.i tends to 1 with time, the acoustic properties of the 
medium are progressively recovered in a slow healing pro- 
cess. The parameter R n , n -i in Fig. 0(b) shows the level 
of resemblance of two successive signals. It is interest- 
ing to see that in this healing process at the end of the 
protocol, two successive signals are practically identical. 
The process of healing is slowly cumulative. This effect 
can be seen as a slow dynamic effect observed for elas- 
tic waves in other materials too, such as sandstones and 
cracked solids [H, Q : the strong wave action modifies 
the acoustic properties of the medium, which arc slowly 
recovered after some time, long compared to the wave 
period. The probing of this slow process is sometimes a 
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powerful indication of damage (associated with the pres- 
ence of internal contacts) in materials (4f| . 

The last parameter plotted in Fig. [TJc) is the charac- 
teristic time T e corresponding to the arrival time of half 
the total energy of the signal. It exhibits a quite clear 
amplitude dependence and has the same qualitative be- 
havior as the parameter R n ^. Until the 30th amplitude 
of the protocol (£33 ~ 0.3 x 10 -5 ), the characteristic time 
T e of the probe wave (even numbers of signals) remains 
unchanged, while the signals sn to S31 exhibit a dimin- 
ishing characteristic time T e . It means that cither the 
wave packet is propagating faster (nonlinear hardening 
of the medium), either the signal energy is preferentially 
attenuated at later times. Due to the strong nonlinear 
attenuation observed in the temporal signals of Fig.[6jb), 
the second process seems more adequate to explain the 
T e dependence on amplitude. This is visible in Fig. [61(b) 
when comparing signals sio and sn or sio and S19. No 
shift in time is detectable while the last part of the signals 
are strongly attenuated relatively to their first part. 

In order to analyze these amplitude dependent effects 
on the signal energy attenuation, it is adequate now to 
focus on the energetical properties of the signals. In 
Fig. [8jd), the experimental amplitude protocol is re- 
called. The total energies of the high-pass filtered sig- 
nals C^^(0) are normalized by a quantity proportional 
to the excitation energy £ 2 both for the low excitation 
amplitudes of the protocol (Fig. [U(a)) and for the higher 
excitation amplitudes (Fig. [5Jb)). The parts comprised 
between the dashed-lincs compare the same range of nor- 
malized energy, from normalized values 0.8 to 1.1 roughly, 
that are the extreme values taken by the signals excited 
with the weak probe amplitude. The lowest value ob- 
served in Fig. [H](a) corresponds to an excitation ampli- 
tude range where the parameter R n> \ begins to decrease 
from 1 to a lower value. It is important to note that for 
the normalized energies associated with the higher exci- 
tation amplitudes, in Fig. [5£b), all the values are below 
the region delimited by the dashed lines, which moans 
that energy absorption is always stronger for e a > £1 
even when the acoustic response of the medium has been 
modified by a strong wave action (even signals S40 to S50 
for instance). The general behavior of the normalized 
energy for increasing excitation amplitudes is decreasing 
and monotonous: the higher the excitation amplitude 
is, the higher the energy absorption is. Its associated 
energy is consequently deviating from the e\ law. The 
nonlinear energy of the LF demodulated part is increas- 
ing monotonously with the excitation amplitude. The 
LF energy exceeds the HF energy of the signal from the 
excitation amplitude £41, as seen in Fig. Etc) when the 
LF energy is above the region delimited by the dashed 
lines. 

In order to analyze more precisely the nonlinear self- 
action process for the HF components and the LF en- 
ergy nonlinear increase process, their normalized energies 
are plotted as a function of the excitation amplitude in 
log-log scale, and compared to power laws in Fig. [5J In 
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FIG. 8: (a) Energy of the high-pass filtered signals C„ }71 (0) 
normalized by a quantity proportional to the excitation en- 
ergy for the low excitation amplitudes of the protocol (d). 
(b) Energy of the high-pass filtered signals for the higher exci- 
tation amplitudes of the protocol, (c) Energy of the low-pass 
filtered signals for the higher excitation amplitudes of the pro- 
tocol, (d) Protocol for the successive excitation amplitudes 
of the experiment. 



Fig. HJa) , the normalized HF energy is plotted as a func- 
tion of the excitation amplitude. This is the same quan- 
tity as that of Fig. [H](b) . For the lowest excitation ampli- 
tudes, the normalized C^n (®) / £ n is decreasing close to 
a power law £~ ' 8 , which means that the energy C^(0) 
is increasing with a £„ dependence. In a perfectly lin- 
ear system, it would be a e\ dependence. Note that the 
overall behavior is close to a £q 25 dependence, which is a 
dramatic deviation from the linear law £ 2 and also a very 
strong self-action compared to usually observed satura- 
tion processes in rocks for instance [461 ] . For higher exci- 
tation amplitudes, a small plateau is observed, denoting 
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a local comportment law close to the linear one e 2 , even if 
the process leading to this locally linear behavior may be 
due to combined effects of different nonlinear processes: 
self-transparency, slow-dynamics, conditioning... For the 
highest excitation amplitudes, the dependence on ampli- 
tude is the strongest, approaching a e°' 5 law for Canity- 
This means that if the excitation energy is multiplied by 
a factor 4 for instance, the received energy is multiplied 
by a factor y/2. This would be an even stronger effect for 
the late part of the signal considering the fact that most 
of the nonlinearly attenuated energy is the one arriving 
at the latest time in the signal, according to Fig. [5] and 
Fig. [Jj Surprisingly, the largest amplitudes of the signal 
are not strongly nonlinearly attenuated but the late con- 
tributions are. This peculiar observation is discussed in 
section Mil 

In Fig.[5]Jb), the normalized LF part energy C% ^(0)/e 2 
is plotted in log-log scale as a function of the excita- 
tion amplitude. A quite clear power law £*' 27 is found 
for this NL process of frequency-down conversion (self- 
demodulation), which corresponds to a LF energy pro- 
portional to £ 3 ' 27 . The classical NL behavior observed in 
homogeneous media with a quadratic nonlincarity is e\ 
for the LF energy [3l|. The observed dependence corre- 
sponds to a Eq 63 power law for the field itself, to be com- 
pared to the e 2 law classically observed. This dependence 
is also closely observed with another processing method 
applied to the LF part: the normalized peak-to-peak 
amplitude of the LF contribution e^ F /e 2 is plotted in 
Fig-OJc) in log-log scale. A power law behavior e~ °' 31 can 
be observed, that is a law e^ 69 for the dynamics of e^ F . If 
compared to the behavior obtained with C^ F (0) oc e 3 ' 27 , 
taking into account that C^ F (0) oc (s!^ H ) 2 , this dynam- 
ics corresponds to a £i- 69x2 = e 3 38 law. This is an in- 
termediate dependence between the classical one e a and 
the Hertzian clapping limit e 3 . The latter dependence 
law has been observed earlier in different experimental 
conditions in Q. Note that for symmetry reasons, the 
Hertz-Mindlin hystcrctic nonlincarity is not able to pre- 
dict the effect of self-demodulation observed here at the 
leading order of the nonlinear approximation, and may 
not contribute [13, S|| . 
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FIG. 9: (a) Energy of the high-pass filtered signals C^^(0) 
normalized by a quantity proportional to the excitation en- 
ergy e\ as a function of the excitation amplitude in log-log 
scale, (b) Energy of the low-pass filtered signals C^„(0) nor- 
malized by a quantity proportional to the excitation energy 
as a function of the excitation amplitude in log-log scale (c) 
Peak-to-peak amplitude of the self-demodulated signal nor- 
malized by in log-log scale. 



III. SUMMARY OF THE RESULTS AND 
DISCUSSIONS 

Different experimental observations have been pre- 
sented for the first time and the main results can be 
summarized as the following: 

• observation of the coexistence of a nonlinear self- 
demodulated LF wave together with a HF coda- 
type wave (Fig. [2]), 

• a drastic increase of the HF coda energy with static 
stress, compared to the dependence of the LF wave 
energy (Fig. HI), 



• a strong shape dependence of the coda on the 
acoustic excitation amplitude (Fig. [5] and Fig. [6]), 
and especially a strong nonlinear attenuation of the 
coda wave at the latest arrival times, 

• almost no variation in the acoustic properties of the 
medium after moderate acoustic excitations up to 
an excitation acoustic strain of e a ~ 0.3 x 10~ 5 , cor- 
responding to signal number 31 (this represents a 
ratio of dynamic over static deformation of s a /£o — 
5 x 1CT 3 ), 

• dynamic modifications of the acoustic properties 
(nonlinear self-action, fast dynamics) even at low 
excitation levels e a ~ 2x 10 -7 corresponding to sig- 
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nal number 11 (this represents a ratio of dynamic 
over static deformation of e a /eo ~ 3 x 10 -4 ), 

• modification of the acoustic properties after large 
excitation levels (slow dynamics and conditioning 
effects), followed by a slow recovery. 

Having in mind the ray acoustic approximation in this 
experiment of multiple scattering of elastic waves, the 
first part of the coda signals can be associated with short 
propagation paths while the last part of the coda sig- 
nals may be associated with longer propagation paths. 
Note that the lengths of the long propagation paths can 
exceed the container size by a factor of 10, which means 
that reflexions on the container sides may occur. In prin- 
ciple, nonlinear attenuation should be more visible on the 
long propagation paths than on the short ones because 
the distance for the accumulation of this nonlinear pro- 
cess is larger [3l|, H?J . This explains the stronger effect 
of nonlinear attenuation on the last part of coda signals, 
that experienced longer propagation paths. The physi- 
cal nature of this nonlinear attenuation process remains 
undetermined. Nonlinear scattering of the elastic waves 
due to opening or closing of the weakest contacts and 
therefore dynamic switching of propagation paths under 
the wave action could explain the nonlinear attenuation. 
This is in agreement with the drastic dependence of the 
HF coda wave energy transmission on the applied static 
stress. Nonlinear absorption of the elastic waves at the 
contacts by stick-slip [42j or thermo-elasticity [43| pro- 
cesses could also contribute. Recent numerical simula- 
tions of the acoustic propagation through unconsolidated 
granular packings [43 . IfjOj could provide an interesting 
insight in the processes involved in the experimentally 
observed nonlinear attenuation. 

For both nonlinear scattering and nonlinear absorp- 
tion, the weakest contacts are supposed to be responsi- 
ble for the main contribution, especially for such small 
acoustic excitation strains as e a ~ 10~ 7 . Coda signals are 
shown to be particularly sensitive to these weak contacts, 
certainly due to the fact that energy is well distributed 
over the entire medium, including weak contacts, after 
few milliseconds (coda signals have durations larger than 
2 ms, which corresponds to propagation paths much 
larger than the emitter - receiver distance). Moreover, 
the wavelength being of the order of the bead diameter 
(A ~ d), small features as small as the beads and their 
contacts are resolved by the acoustic wave. This is not 
the case for long wavelength propagation (A ^> d), which 
is known to be less sensitive to the weak contacts in the 
medium and more influenced by the average properties 
of the medium 0, . 

The modification of the acoustic properties after the 
application of a strong excitation, is commonly re- 
lated to the effect of conditioning and slow dynamics 
[HI, IH, El, HH ■ This addresses the problem of the acous- 



tic wave influence on the granular medium state and 
memory. Under a given threshold, the medium is appar- 
ently not modified (under excitation amplitude number 
31, e a ~ 0.3 10~ 5 which corresponds to a ratio of dynamic 
over static deformation of e a /eo ~ 5 x 10 -3 ), while over 
some excitation amplitude, the medium exhibits mem- 
ory of the past acoustic excitation. Coda signal analysis 
is shown to be a powerful tool to monitor slight mod- 
ifications in the acoustic response of an unconsolidated 
granular structure as a function of time, the resolution 
(including some time averaging) being of the order of the 
second. 



IV. CONCLUSIONS 

In this work, experimental results on the nonlinear 
transport of short wavelength acoustic waves (of the or- 
der of the bead diameter) are reported. It is shown 
that the amplitude dependence of the codas is strong, 
and contains some information on the strength of the 
propagation paths in the medium. An amplitude depen- 
dent attenuation is visible on the characteristic time T e , 
which diminishes with the excitation amplitude: most of 
the nonlinear self-action effects take place at the latest 
arrival times of the coda signals. This shows that the 
field which encountered the larger number of scattering 
events, or which does not follow the strongest force chains 
is more influenced by the effects of nonlinear dissipation 
(stick-slip, sliding, or thermo-elastic losses) and nonlin- 
ear scattering (closing or opening of the weak contacts 
for instance). 

Among the possible applications of such experimental 
results is the non destructive testing (NDT) of granular 
structures with nonlinear acoustic methods, in particular 
with multiple scattered waves (coda signals), in contrast 
to the widely used and recent nonlinear methods employ- 
ing coherent wave interactions [H, . For this purpose, 
it would be of interest to make use in the future of the 
coda wave interferometry technique [H, [5f| [56| . 

It has been shown that the acoustic properties of the 
medium can be modified temporary by a strong acous- 
tic excitation. In the context of NDT, this work show 
the range of acoustic amplitudes to use (or the time nec- 
essary to recover the initial state) in order to stay in a 
regime where the acoustic waves do not modify the clastic 
response of the medium. 
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